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ABSTRACT 

We have used the Very Large Array(VLA) in C configuration to carry out a 
sensitive 20cm radio survey of regions of sky that have been surveyed in the Far Infra- 
Red over the wavelength range 5-200 microns with ISO as part of the European Large 
Area ISO Survey (ELAIS). As usual in surveys based on a relatively small number of 
overlapping VLA pointings the flux limit varies over the area surveyed. The survey 
has a flux limit that varies from a Scr limit of 0.135mJy over an area of 0.12deg^ to a 
5cr limit of 1.15mJy or better over the whole region covered of 4.22 deg^. In this paper 
we present the radio catalogue of 867 sources. These regions of sky have previously 
been surveyed to shallow flux limits at 20cm with the VLA as part of the VLA D 
configuration NVSS(FWHM=45 arcsec) and VLA B configuration FIRST(FWHM=5 
arcsec) surveys. Our whole survey has a nominal 5 sigma flux limit a factor of 2 below 
that of the NVSS; 3.4 deg^ of the survey reaches the nominal flux limit of the FIRST 
survey and 1.5 deg^ reaches to 0.25 mJy, a factor of 3 below the nominal FIRST survey 
limit. In addition our survey is at resolution intermediate between the two surveys and 
thus is well suited for a comparison of the reliability and resolution dependent surface 
brightness effects that affect interferometric radio surveys. We have carried out a a 
detailed comparison of the reliability of our own survey and these two independent 
surveys in order to assess the reliability and completeness of each survey. 

Key words: radio continuum: general - galaxies 



1 INTRODUCTION 

Tlie Infrared Space Observatory (ISO, Kessler et al. 1996), 
launched in November 1995 was the successor of the Infrared 
Astronomical Satellite (IRAS) and provided unparallel sen- 
sitivity in mid to far infrared wavelengths (i.e. 5-200 /im). 
The European Large- Area ISO Survey (ELAIS, Oliver et al. 
1997, Oliver et al. 1998 in preparation) is a project that used 



ISO to carry out a deep wide angle survey at wavelengths of 
6.7, 15, 90 and 175 /im. The 6.7 and 15 /im surveys were car- 
ried out with the ISO-CAM camera (Cesarsky et al. 1996) 
with the aim to reach a 5(t sensitivity of ~2mJy at 15mi- 
crons. The 90^m and 175fim surveys used the ISO-PHOT 
camera (Lemke et al. 1994) with the aim to reach a 5a sen- 
sitivity of ~25mJy. At these limits, we expect ISO to be 
confusion limited at 90 fim and 175 jj,m by galaxies and 
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galactic cirrus emission and hence this survey should be the 
deepest FIR survey possible with the satellite. 

The area covered in the ELAIS survey is ~ 13 square 
degrees at 15 and 90 microns, ~ 7 square degrees at 6.7 
microns and ~ 3 square degrees at 175microns. 

The ELAIS survey is ~50 times deeper at 5-20/im 
than IRAS. Thus our survey will allow us to explore IRAS- 
like populations to higher redshift and possibly unveil new 
classes of objects or unexpected phenomena. We expect to 
detect thousands of galaxies, many of which will be at high 
redshifts and undergoing vigorous star formation. The ex- 
pected large number of high-z IR galaxies should provide 
vital information about the star formation rate out to z=l 
and possibly earlier. 

The spatial resolution of ISO will be insufficient to prop- 
erly identify optically faint objects. At 15 microns, the sur- 
vey resolution is ~ 10 arcsec and at 90 microns it will be 
about one arc minute. Complementary radio data will play 
a crucial role in identifying many of the most interesting ob- 
jects, as they did in the early days of X-ray astronomy (e.g. 
Cyg X-1) and in more recent times for IRAS (e.g. IRAS 
F10214+4714 (Rowan- Robinson et al. 1991). 

In this paper we report the description of the radio 
observations obtained in the three ISO-ELAIS survey re- 
gions in the northern celestial hemisphere (Nl_1610+5430, 
N2_1636+4115 and N3_1429+3306). The observations are 
made with the Very Large Array (VLA) radio telescope at 
1.4GHz (20cm) in the VLA C-configuration (maximum base- 
line 11km) with a resolution (FWHM) of ~ 15 arcsec. The 
aim of these VLA observation was to obtain an uniform 
covering of the ELAIS regions with a rms noise limit of ~50 
/xJy. These VLA observations will be essential in the optical 
identification phase of the ELAIS sources and in assessing 
the reliability of the ELAIS source lists. 

Moreover, with a radio survey it will be possible to in- 
vestigate the radio/far-infrared correlation in star forming 
galaxies to flux densities deeper than those reached by IRAS. 
Helou, Sofier & Rowan-Robinson (1985) noted a strong cor- 
relation between radio and far infrared flux for star form- 
ing galaxies, valid over a very wide range of infrared lumi- 
nosities, and this has been confirmed in many other studies 
(e.g. Wunderlich, Klein & Wielebinski 1987; Condon, An- 
derson & Helou 1991). The radio emission is interpreted as 
a synchrotron radiation from relativistic electron which have 
leaked out of supernova remnants. It is expected that this 
correlation should extend below the IRAS flux level since 
the majority of the sub-mjy radio sources have been iden- 
tified with faint blue galaxies whith spectra similar to those 
of star forming objects (Benn et al. 1993). 

In addition, combining deep radio and optical data with 
the ISO survey fluxes will provide information on the trivari- 
ate IR-radio-optical luminosity function and its evolution 
and the contribution of starburst galaxies to the sub-mJy 
radio source counts. The ratio of the FIR emission and ra- 
dio emission will also allow is to investigate the physical 
origin and spatial distribution of the energy sources in the 
detected objects in the same way that VLA maps have been 
central to our understanding of the origin of IRAS sources. 

Finally, this survey, due to its depth and extension, is 
very important also as radio survey in its own right. In fact, 
the selected sample is large and deep enough to constitute 
a statistically significant sample of sub-mJy radio sources. 



whose nature and characteristics are still a major topic in ob- 
servational cosmology (see Windhorst, Mathis & Neuschae- 
fer 1990, Fomalont et al. 1991, Rowan- Robinson et al. 1993, 
Gruppioni et al. 1997) 



2 RADIO OBSERVATIONS 

2.1 Choice of observing frequency and VLA 
configuration 

The VLA C-configuration and the observing frequency of 
1.4 GHz give the optimum resolution to acquire the kind of 
radio data that we need. Whilst less prone to surface bright- 
ness effects , the VLA D configuration is confusion-limited 
at the fluxes we wish to attain (the 5 a confusion limit in 
D conflguration is 0.4 mjy/beam). With the C conflgura- 
tion and a frequency of 1.4 GHz the synthesized beam size 
(Full Width at Half Power, FWHP) is ~15 arcsec. The well- 
defined synthesized beam of the VLA should enable us to 
pinpoint optical identifications to 1 arcsec, except for the 
asymmetric multi-components sources. The frequency of 1.4 
GHz was chosen because at this frequency the FWHP of the 
VLA primary beam is 31 arcmin. This allow us to cover the 
ELAIS field with a relative small number of pointing cen- 
ters. In fact it is possible to obtain a mosaic map with nearly 
uniform sensitivity if the separation is 31 / \/2 ~ 22 arcmin. 
Moreover, at 1.4 GHz there will be contributions from both 
the steep and fiat spectrum population of radio sources. 

Our observations are made in spectral line mode us- 
ing two different IF channels centered at 1.3649 GHz (IFl) 
and 1.4352 GHz(IF2). Each IF channel has a bandwidth of 
18.75 MHz, which is subdivided into 7 spectral line chan- 
nels evenly spaced in frequency across the bandwidth of the 
input IF channel. Therefore, we have a total of 14 spectral 
line channels (7 for each IF) with a bandwidth of 2.68 MHz 
each. We decided to use the line mode in order to facilitate 
wide field mapping and avoid the effect of interference. The 
total bandwidth available is thus 37.5 MHz, narrower than 
the 50MHz used in continuum mode. A narrower bandwidth 
means a worse sensitivity. In our case we lose about ~25% 
in sensitivity or 15% in signal to noise near the pointing 
centers because the total bandwidth is 75% of that in the 
continuum mode. However, we considerably reduce the chro- 
matic aberration (bandwidth smearing), which reduces the 
area covered by each pointing in continuum mode. In addi- 
tion, the line mode is less susceptible to narrow interference 
noise spikes, since one only needs to excise the channel that 
is affected rather than loosing the whole IF band, which is 
the case in continuum mode. Since our aim is to obtain an 
uniform sensitivity over the whole ELAIS fields and not to 
obtain a single-field deep survey, we have opted to use the 
line mode. Moreover, the reduced bandwidth smearing in 
line mode will give us more accurate angular sizes for our 
sources. 



2.2 Observations 

The VLA observation of the ELAIS regions were carried out 
in Aprfl 96 (10 hours) and in July 1997 (24 hours). A total of 
20 pointings were made in 20cm spectral line mode. Seven of 
these pointings are in the ELAIS field Nl, eight in the field 
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Figure 1. The sky position and orientation of the ISO survey 
regions. Each outer rectangle is 120' X 80'. The triangles show 
the ISO CAM pointings, while the asterisk show the ISO PHOT 
pointings. The circles show the regions we mapped with the VLA. 
The circles are drawn with radii R= 20', where the VLA power 
sensitivity is ~30% of the central value. The thick circle in N2 is 
the deep pointing N2 VD (see text for more details). 



N2 and five in N3. Tlie integration time of each pointing 
was ~ 1 liour. This allowed us to obtain, in each field, a 
root mean square (rms) noise of ~ 0.05 mjy (a 5 ct limit of 
0.25 mJy). Moreover, during the July 1997 run, we observed 
two regions (one in Nl and one in N2) with an integration 
time of ~ 3 hours each. In particular, in Nl we re-observed 
(for 3 hours) the pointing VI (already observed for one 
hour in April 1996) while in N2 the deep pointing (N2 VD) 
was shifted of 10 arc min from the map center to avoid the 
presence of a strong (>100 mJy) radio source. In the two 
deep pointings Nl VI and N2 VD we reached an rms noise 
of ~ 0.026 mJy (a 5 a limit of 0.13 mJy). Finally, in order 
to study the degradation of the image quality as the off- 
axis of the sources increase, a strong calibrator (the source 
3C286) was observed at position offset of 5,10,15,20,25,30 
and 35 arcmin in two orthogonal directions (North-South 
and E ast-West). The result of this test is reported in section 



5.1.2 



In Table |l| we report the position of the center of each 
observation, while Figure ^ shows the sky position and ori- 
entation of the ISO survey regions. The circles show the 
VLA regions mapped. The circles are drawn with a radius 
(R~ 20') where the VLA power sensitivity is ~30% of the 
central value. 



3 DATA REDUCTION 

All the data were analyzed with the NRAO AIPS reduc- 
tion package. The data were calibrated using 3C286 as pri- 
mary flux density calibrator (assuming a flux of 15.04 Jy 
at 1.3649 GHz (IFl) and 14.70 Jy at 1.4352 GHz (IF2)) 
and the sources 1549-1-506 and 1635-1-381 as secondary cal- 
ibrators. The task UVFLAG was used to "flag" (delete) the 
corrupted data (e.g. bad integration, non operating anten- 
nas, high amplitudes due to interferences etc) while the tasks 
VLACALIB and GETJY were used to calibrate the data and to 
determine the source flux densities. Finally each observation 
was cleaned using the task IMAGR. 



3.1 Root mean square (rms) map noise of the 
single pointings 

The integration time of each observation after deletion of 
the corrupted data are reported in Table ^. The rms noise 
of each pointing was estimated using the amplitude distri- 
bution of the pixel values in the cleaned map before cor- 
recting for primary-beam attenuation. In Figure ^ we show 
this distribution for the pointing N2 VI. This distribution is 
the sum of a Gaussian noise core with a positive-going tail 
caused by discrete sources. Note that the y-axis in Figure ^ 
is logarithmic. 

Thus the rms of the noise distribution alone should be 
nearly equal to the rms of that distribution obtainable by 
reflecting the negative flux portion of the observed ampli- 
tude distribution about flux=0. In column VII of Table ^ 
we report the theoretical rms noise level (at la) as com- 
puted directly from the observing parameters (integration 
time, number of antennas, observing frequency, bandwidth 
and number of IF channels) , while in column VIII we report 
the Ict rms noise level obtained by fitting the noise distri- 
bution of each pointing. As shown in Table ^ there is a very 
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Table 1. VLA Observations of ELAIS regions 



Region Pointing Observing RA DEC Integration Theoretical Observed 

Data (J2000) (J2000) Time (min) rms (mjy) rms (mjy) 



Nl 


VI 


April 


96 


16 


10 


00.00 


+54 


30 


00.00 


58.4 


0.052 


0.049 


Nl 


VI 


July 


97 


16 


10 


00.00 


+54 


30 


00.00 


175.4 


0.030 


0.030 


Nl 


V2 


July 


97 


16 


09 


25.51 


+54 


06 


32.42 


57.8 


0.052 


0.050 


Nl 


V3 


July 


97 


16 


07 


34.62 


+54 


25 


05.31 


58.5 


0.052 


0.051 


IN JL 


V 4 


July 


y ( 


ID 


Uo 


no 

uy .DO 


+04 




Uo. 14 


oy.o 




U.UOz 


Nl 


V5 


July 


97 


16 


10 


37.76 


+54 


54 


36.12 


59.0 


0.052 


0.049 


Nl 


V6 


July 


97 


16 


12 


28.43 


+54 


35 


55.56 


59.0 


0.052 


0.051 


Nl 


V7 


July 


97 


16 


11 


51.23 


+54 


11 


54.45 


57.8 


0.052 


0.052 


N2 


VI 


April 


96 


16 


36 


00.00 


+41 


06 


00.00 


65.3 


0.050 


0.050 


N2 


V2 


April 


96 


16 


34 


00.00 


+41 


06 


00.00 


56.8 


0.053 


0.052 


N2 


V3 


April 


96 


16 


38 


00.00 


+41 


06 


00.00 


59.3 


0.052 


0.050 


N2 


V4 


July 


97 


16 


37 


00.00 


+40 


46 


00.00 


59.2 


0.052 


0.052 


N2 


V5 


July 


97 


16 


37 


00.00 


+41 


26 


00.00 


59.0 


0.052 


0.052 


N2 


V6 


July 


97 


16 


35 


00.00 


+40 


46 


00.00 


59.5 


0.052 


0.053 


N2 


V7 


July 


97 


16 


35 


00.00 


+41 


26 


00.00 


59.5 


0.052 


0.051 


N2 


VD 


July 


97 


16 


35 


00.00 


+41 


06 


00.00 


178.4 


0.030 


0.050 


N3 


VI 


April 


96 


14 


28 


30.00 


+33 


20 


00.00 


56.5 


0.053 


0.053 


N3 


V2 


April 


96 


14 


26 


30.00 


+33 


20 


00.00 


58.7 


0.052 


0.053 


N3 


V3 


July 


97 


14 


27 


30.00 


+32 


58 


30.00 


58.9 


0.052 


0.052 


N3 


V4 


July 


97 


14 


30 


30.00 


+33 


20 


00.00 


59.5 


0.052 


0.050 


N3 


V5 


July 


97 


14 


29 


30.00 


+32 


58 


30.00 


58.0 


0.052 


0.051 




16 12 11 10 09 

RIGHT ASCENSION (J2000) 

Figure 2. Logarithmic histogram of VLA map peak flux densities 
in the pointing N2 VI 



good agreement between the observed and the theoretical 
rms noise levels. 



Figure 3. Central 26x26 arcmin^ = 520x520 pixels of the mosaic 
map Nl. The contour level are 0.25,0.50,0.75,1.00, 1.25 ... mjy. 
The noise at the corners of the image (see for example the low 
right corner) is due to the increase of the rms noise as the off-axis 
value increase 



© 1998 RAS, MNRAS 000, 



A Deep VLA survey at 20cm of the ISO ELAIS survey regions 5 



3.2 Mosaic maps 

Using the AIPS task LGEDM, HGEOM and LTESS we have com- 
bined all the observations and we have created a mosaic map 
for each field (Nl, N2 and N3). A special procedure has been 
adopted for N2. The presence of a deep observation not in 
the center of the radio map (see above and Figure ^ makes 
the noise of this map strongly irregular. To simplify the ex- 
traction of the sources we have created two different mosaic 
maps. In the first map we have combined all the observation 
in N2 excluding the deep pointing, while in the second one 
we have combined the deep pointing N2 VD with all the sur- 
rounding pointings (N2 VI, N2 V2, N2 V6 and N2 V7, see 
Figure hi). In this way we have obtained two mosaic maps 
with a regular noise: lower in the map center and higher in 
the outer regions. 

An example of our map is given in Figure ^ where we 
show the contour map of the central region (26x26 arcmin^) 
of the mosaic map Nl. 



3.3 Noise of the mosaic maps 

We analyzed the noise properties in each mosaic maps (Nl, 
N2, N2 Deep and N3). As expected, Nl, N2 and N2 Deep 
have a regular noise distribution: a circular central region 
with a flat noise surrounded by concentric annular region, 
where the noise increases for increasing distance from the 
center (the off-axis value). No structures or irregularities 
were found in the rms maps. In Figure ^we plot the varia- 
tion of the rms as function of the off-axis value for Nl, N2 
and N2 Deep. Due to the different number of pointings in 
N3 (five pointings instead of seven, see Figure |l|) the mo- 
saic map in this field has a different shape and the noise 
distribution in the map has a semi-circular shape instead of 
a circular one. In Figure |^,e,f we plot the variation of the 
rms value in three different slices. In Figure kli we plot the 
rms in the North-South direction, in Figure He the rms in a 
slice that connects the center of the three northern pointings 
(N3 VI, N3 V2 and N3 V4), while in Figure |f we report 
the rms value in a slice that connects the center of the two 
souther pointings (N3 V3 and N3 V5). As for Nl and N2, 
no structures or irregularities were found in the rms map of 
N3. 



4 SOURCE CATALOGUE 

4.1 Regions for the source extraction 

Using Figure ^ we determined the regions with constant rms 
noise to be used for the source extraction. Due to the circu- 
lar shape of the rms map, we searched for radio sources in 
a circular central region plus concentric annular regions to 
an off-axis value of 42 arcmin. Outside this area, the beam 
attenuation significantly increases the limiting fiux (see Fig- 
ure ^ . In each region the rms used for source extraction is 
the highest rms in the region. In this way we are sure that 
all the sources extracted are above the fixed threshold in 
the region (for example 5 times the rms value). The sizes of 
these regions are reported in Table ^ and in Figure ^ 

Using the data reported in Table ^ we obtained the 
solid angle Q{Sp) as function of the peak fiux density Sp 



Table 2. The regions used for the sources extraction 



Region 




Rout 


Size 


Total 


rms 


5 (T limit 




(') 


(') 


deg^ 


dcg2 


mJy 


mJy 


Nl 














1 





8 


0559 


0559 


0.027 


0.135 


2 




14 


0.1151 


0.1710 


0.034 


0.170 


3 


14 


22 


0.2513 


0.4223 


0.038 


0.190 


4 


22 


26 


0.1677 


0.5900 


0.050 


0.250 


5 


26 


30 


0.1955 


0.7855 


0.070 


0.350 


6 


30 


34 


0.2234 


1.0089 


0.100 


0.500 


■7 
1 


o4 


QQ 
OO 


V.ZOi-O 


1 ofino 
i.zDUz 


U. iOU 


U. / OU 


8 


38 


42 


0.2 f\)o 


1.5395 


0.230 


1.150 


N2 














1 





26 


5000 


5900 


0.050 


0.250 


2 


26 


30 


0.1955 


0.7855 


0.070 


0.350 


3 


30 


34 


0.2234 


1.0089 


0.100 


0.500 


4 


34 


38 


0.2513 


1.2602 


0.150 


0.750 


O 


Oo 


42 


n 97Q'\ 
u.z ( yo 


1 5*^95 


u. zou 


1 1 5n 


N2 D 














1 





8 


0.0559 


0.0559 


0.027 


0.135 


2 


8 


11 


0.0497 


0.1056 


0.030 


0.150 


3 


11 


14 


0.0655 


0.1711 


0.038 


0.190 


N3 














1 





26 


0.2950 


0.2950 


0.050 


0.250 


2 


26 


30 


0.1645 


0.4595 


0.070 


0.350 


3 


30 


34 


0.1961 


0.6556 


0.100 


0.500 


4 


34 


38 


0.2279 


0.8835 


0.150 


0.750 


5 


38 


42 


0.2597 


1.1432 


0.230 


1.150 



covered by our survey. This effective survey area is shown in 
Figure 0. 



4.2 Source detections 

Within each region we searched for radio sources up to a 
peak flux density > 5 times the rms value of the region. 
The sources were extracted using the task SAD (Search And 
Destroy) which attempts to find all the sources whose peaks 
are brighter than a given level. For each selected source the 
flux, the position and the size are estimated using a least 
square Gaussian fit. However, for faint sources a Gaussian 
fit may be unreliable (see Condon 1997, for an extensive 
discussion about errors in Gaussian fits). For this reason 
we decided first to run the task SAD with a flux limit of 
3 times the local rms value (see Table |^. Subsequently, 
we derived the peak flux of faint sources (those detected 
with 3(T < Speak < 7a) using a second degree interpolation 
(task MAXFIT). Only the sources with a MAXFIT peak flux 
density > 5 a were included in the sample. For these faint 
sources the total flux density was obtained using the task 
IMEAN, which integrates the map value in a specific rectangle, 
while for all the other parameters (major axis, minor axis 
and position angle) we used the values obtained with the 
Gaussian fit. For irregular resolved sources the total flux 
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Figure 4. The rms noise value as function of the off-axis value. The solid line represent an rms value of 0.05 mjy. See text for more 
details. 



density was calculated using the task TVSTAT which allows 
us to use an irregular area to integrate the map value. 



4.3 Multiple sources 

Detected sources separated by less than two times the value 
of our synthesized beam size (i.e. 25 arcsec) and with a flux 
ratio lower than 2 have been considered as an unique source. 
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Figure 5. The region used for the source extraction. The Unes with the same thickness enclose regions with the same rms noise value 
(see Table H 



We adopted this criterion because the component flux den- 
sity ratio of physically doubles is usually small ( S 2) while 
the projection pairs can have arbitrarily large flux density 
ratio (Condon, Condon and Hazard, 1982). 



5 THE SOURCE CATALOGUE 

Considering all the available observations we detected a total 
of 867 sources at > 5 cr level (44 of which have multiple 
components) over a total area of 4.222 deg^. The catalogue 
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321 
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Figure 7. Peak flux distributions (upper) and total to peak flux 
ratio as function of peak flux (lower) for all the radio sources 
detected in Nl, N2 and N3 



with all the 867 sources (921 components) reports the name 
of the source, the peak fiux density Sp in mjy, the total 
flux density S/ in mJy, the RA and DEC (J2000), the fuU 
width half maximum (FWHM) of the major and minor axes 
Om and 6m (in arcsec), the positional angle PA of the major 
axis (in degrees) and the off-axis values in the VLA map 
(in arcmin). The different components of multiple sources 
are labeled "A", "B", etc., followed by a line labeled "T" in 
which flux and position for the total sources are given. For 
these total sources the position have been computed as the 
flux- weighted average position for all the components. Table 
3 shows the first page of the catalogue as an example. 

In Table |^ we report the number of radio sources de- 
tected in Nl, N2 and N3, while in Figure ^ we show the 
distribution of the peak fiux density and the total to peak 
flux ratio as a function of peak flux for all the 867 sources. 
Contour maps of all the 44 double or multiple sources are 
shown in Figure ^ ^ and |lo| 

5.1 Sources parameters and their uncertainties 

5.1.1 Positions and angular sizes 

Following Condon 1997, the error on the sources parameters 
reported in the catalogue have been calculated using : 



SI S| 91, el, 2 91,61, ) p2 ^ ^ 
where p is the signal to noise ratio given by 

_ 71" 6M6mSp , , 

^ " 8ln2 h^alap ^ ' 

Equations 1 and 2 are the master equations for estimating 
variances of the parameters derived from a two dimensional 
Gaussian flt on an image with a noise variance a^ap and 
pixel size h. In our maps h=3 arcsec, while column 5 of 
Table ^ reports the value of amap as function of the off-axis 
value. The rms position error are given by (Condon et al. 
1998): 

(J^ = £q 4- (Jej^j sin^ (PA) + al^ cos^ {PA) (3) 

aj = el + ol^, cos'' {PA) + al^ sir? {PA) (4) 

The mean image off-set < Aa >, < /S.S > and rms calibra- 
tion uncertainties and el axe best determined by compar- 
ison with accurate position of sources strong enough that 
the noise plus confusion terms are much smaller than the 
calibration terms. We used the 37 single compact sources 
stronger than 10 mJy found in the FIRST (Faint Images 
of the Radio Sky at Twenty cm) radio survey (Becker et al. 
1997). Their off-sets < Aa > and < A5 > are shown in Fig- 
ure hll The mean off-sets of our map are < Aa —0.15" 
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Table 3: The VLA Radio Catalog of the ISO Regions Nl N2 and N3 



Name Sp asp S/ as, R.A. (J2000) DEC(J200) as 9m 6^™ PA ae,,, (t6I™ apA Off axis 

(mjy) (mjy) (mjy) (mjy) (//) (//) (//) (//) deg (//) (//) deg arcmin 



ELAISR142529+331644 


1.145 


0, 


.057 


1.272 


0, 


.064 


14 


25 


29, 


.3272 


+33 


16 


44, 


.204 


0, 


.837 


0, 


.875 


15. 


.959 


13, 


.577 


-2, 


.2 


0, 


.801 


0, 


.681 


0, 


.2 


37, 


.9 


ELAISR142529+332324 


0.829 


0, 


.056 


0.955 


0, 


.065 


14 


25 


29, 


.8750 


+33 


23 


24, 


.705 


1, 


.279 


0, 


.968 


17, 


.910 


12, 


.515 


70, 


.1 


1, 


.219 


0, 


.852 


0, 


.1 


37, 


.8 


ELAISR142532+330126 


7.325 


0, 


.089 


8.049 


0, 


.097 


14 


25 


32 


.7110 


+33 


01 


25, 


.991 


0, 


.520 


0, 


.390 


15, 


.582 


13, 


.752 


108, 


.7 


0, 


.188 


0, 


.166 


0, 


.1 


41, 


.5 


ELAISR142536+331600 


1.855 


0, 


.054 


2.309 


0. 


.068 


14 


25 


36 


.5493 


+33 


16 


00, 


.729 


0, 


.672 


0, 


.571 


16, 


.556 


14, 


.659 


-52, 


.1 


0, 


.484 


0, 


.429 


0, 


.2 


36 


.5 


ELAISR142544+333704 


1.399 


0, 


.095 


1.317 


0, 


.090 


14 


25 


44, 


.7981 


+33 


37 


04, 


.548 


1, 


.211 


0, 


.858 


16, 


.427 


11, 


.240 


78, 


.7 


1, 


.121 


0, 


.767 


0, 


.1 


38, 


.5 


ELAISR142552+331823 


0.726 


0, 


.043 


0.640 


0, 


.038 


14 


25 


52, 


.0659 


+33 


18 


23 


.576 


1, 


.002 


0, 


.773 


14, 


.816 


11, 


.592 


84, 


.5 


0, 


.878 


0, 


.687 


0, 


.2 


33 


.0 


ELAISR142606+333038 


1.959 


0, 


.045 


3.490 


0, 


.081 


14 


26 


06, 


.1840 


+33 


30 


37 


.886 


0, 


.684 


0, 


.527 


21, 


.413 


16, 


.224 


-69. 


.3 


0, 


.496 


0, 


.376 


0, 


.1 


31 


.8 


ELAISR142610+330008 


1.212 


0, 


.053 


1.610 


0, 


.070 


14 


26 


10, 


.9163 


+33 


00 


08, 


.397 


0, 


.884 


0, 


.747 


17, 


.058 


15, 


.186 


84. 


.6 


0, 


.739 


0, 


.658 


0, 


.3 


35, 


.2 


ELAISR142612+333113 


4.870 


0, 


.053 


6.300 


0. 


.069 


14 


26 


12, 


.1433 


+33 


31 


13 


.778 


0, 


.521 


0, 


.387 


18, 


.000 


14, 


.015 


61, 


.6 


0, 


.197 


0, 


.153 


0, 


.4 


30, 


.9 


ELAISR142613+325535 


1.063 


0, 


.053 


1.402 


0, 


.070 


14 


26 


13 


.6314 


+32 


55 


35, 


.104 


0, 


.900 


0, 


.908 


17, 


.157 


14. 


.990 


-21, 


.1 


0, 


.851 


0, 


.743 


0, 


.3 


37, 


.6 


ELAISR142616+330922 


1.948 


0, 


.027 


2.096 


0, 


.029 


14 


26 


16, 


.7592 


+33 


09 


22, 


.764 


0, 


.528 


0, 


.404 


15, 


.447 


13. 


.583 


50, 


.6 


0, 


.216 


0, 


.190 


0, 


.1 


29, 


.8 


ELAISR142619+332006 


0.806 


0, 


.016 


2.677 


0, 


.052 


14 


26 


19, 


.5985 


+33 


20 


05, 


.840 


0, 


.919 


0, 


.499 


41, 


.726 


15, 


.521 


75, 


.3 


0, 


.803 


0, 


.299 


0, 


.5 


27, 


.2 


ELAISR142624+332015 


1.393 


0, 


.024 


1.941 


0. 


.033 


14 


26 


24, 


.8441 


+33 


20 


15, 


.622 


0, 


.567 


0, 


.445 


17, 


.326 


15, 


.686 


-66. 


.4 


0, 


.298 


0. 


.270 


0, 


.1 


26, 


.1 


ELAISR142626+332147 


0.484 


0, 


.024 


0.696 


0, 


.034 


14 


26 


26, 


.2085 


+33 


21 


46, 


.988 


0, 


.893 


0, 


.985 


19, 


.916 


14, 


.078 


26, 


.3 


0, 


.969 


0, 


.685 


0, 


.1 


25, 


.9 


ELAISR142629+332005 


1.792 


0, 


.018 


2.203 


0, 


.022 


14 


26 


29, 


.0804 


+33 


20 


05, 


.221 


0, 


.514 


0, 


.381 


16, 


.569 


14, 


.468 


85, 


.3 


0, 


.168 


0, 


.147 


0, 


.1 


25, 


.3 


ELAISR142633+333211 


5.102 


0, 


.077 


7.516 


0. 


.113 


14 


26 


33 


.5502 


+33 


32 


11, 


.764 


0, 


.548 


0, 


.433 


16, 


.952 


16, 


.945 


103, 


.7 


0, 


.254 


0. 


.254 


1, 


.5 


27, 


.2 


ELAISR142639+331839 


3.387 


0, 


.019 


3.876 


0, 


.022 


14 


26 


39 


.0383 


+33 


18 


39 


.705 


0, 


.494 


0, 


.360 


15, 


.597 


14, 


.306 


-87, 


.3 


0, 


.087 


0, 


.080 


0, 


.5 


23, 


.2 


ELAISR142640+331505 


0.364 


0, 


.019 


0.407 


0, 


.021 


14 


26 


40, 


.2382 


+33 


15 


05, 


.805 


1, 


.119 


0, 


.693 


19, 


.266 


11, 


.314 


85, 


.1 


1, 


.010 


0, 


.593 


0, 


.1 


23, 


.5 


ELAISR142647+331921 


0.424 


0, 


.018 


0.541 


0. 


.023 


14 


26 


47, 


.3297 


+33 


19 


21, 


.184 


0, 


.802 


0, 


.778 


16, 


.609 


14, 


.987 


18. 


.0 


0, 


.700 


0. 


.631 


0, 


,3 


21, 


.5 


ELAISR142649+325410 


1.178 


0, 


.038 


1.346 


0, 


.043 


14 


26 


49, 


.2893 


+32 


54 


10, 


.400 


0, 


.707 


0, 


.568 


16, 


.154 


13, 


.793 


-75, 


.7 


0, 


.518 


0, 


.442 


0, 


.1 


33 


.3 


ELAISR142650+332940A 


9.921 


0, 


.061 


9.712 


0, 


.060 


14 


26 


50, 


.6308 


+33 


29 


42, 


.645 


0, 


.494 


0, 


.360 


14, 


.725 


12, 


.963 


77, 


.4 


0, 


.091 


0, 


.080 


0, 


.5 


22, 


.9 


ELAISR142650+332940B 


1.430 


0, 


.060 


1.453 


0. 


.061 


14 


26 


48, 


.9130 


+33 


29 


26, 


.737 


0, 


.746 


0, 


.712 


14, 


.725 


13 


.456 


3. 


.3 


0, 


.619 


0. 


.565 


0, 


.3 


23, 


.1 


ELAISR142650+332940T 


9.921 


0, 


.009 


11.045 


0. 


.090 


14 


26 


50, 


.4144 


+33 


29 


40, 


.641 


0, 


.486 


0, 


.351 


























23, 


.0 


ELAISR142652+330051 


1.096 


0, 


.027 


1.231 


0, 


.030 


14 


26 


52, 


.2219 


+33 


00 


51, 


.416 


0, 


.602 


0, 


.508 


15, 


.854 


13 


.814 


-39, 


.0 


0, 


.386 


0, 


.336 


0, 


.1 


28, 


.0 


ELAISR142654+332851 


0.603 


0, 


.055 


0.735 


0. 


.067 


14 


26 


54, 


.3361 


+33 


28 


51, 


.537 


1, 


.917 


1, 


.122 


20, 


.432 


11, 


.622 


85. 


.1 


1, 


.859 


1. 


.057 


0, 


.1 


21, 


.8 


ELAISR142657+331815 


0.400 


0, 


.015 


0.564 


0, 


.022 


14 


26 


57, 


.6414 


+33 


18 


15, 


.688 


0, 


.900 


0, 


.781 


22, 


.214 


15, 


.329 


-51, 


.7 


0, 


.848 


0, 


.585 


0, 


.1 


19, 


.4 


ELAISR142658+333312 


2.516 


0, 


.078 


3.534 


0, 


.110 


14 


26 


58, 


.8540 


+33 


33 


12, 


.066 


0, 


.710 


0, 


.621 


16, 


.659 


16, 


.443 


-72, 


.6 


0, 


.519 


0, 


.512 


1, 


.7 


23, 


.2 


ELAISR142701+325935 


0.964 


0, 


.027 


1.069 


0. 


.030 


14 


27 


01, 


.3063 


+32 


59 


35 


.389 


0, 


.684 


0, 


.494 


17, 


.327 


12, 


.477 


86. 


.9 


0, 


.482 


0. 


.347 


0, 


.1 


27, 


.6 


ELAISR142701+330049 


0.533 


0, 


.024 


0.650 


0, 


.029 


14 


27 


02, 


.0018 


+33 


00 


48, 


.872 


0, 


.893 


0, 


.817 


17, 


.448 


15, 


.461 


48, 


.4 


0, 


.787 


0, 


.697 


0, 


.3 


26, 


.6 


ELAISR142701+330525 


19.359 


0, 


.020 


19.645 


0, 


.020 


14 


27 


01, 


.6296 


+33 


05 


25, 


.073 


0, 


.486 


0, 


.351 


14, 


.514 


13 


.634 


-0, 


.6 


0, 


.015 


0, 


.014 


0, 


.0 


23, 


.5 


ELAISR142701+331222A 


3.712 


0, 


.018 


4.667 


0. 


.023 


14 


27 


02, 


.1196 


+33 


12 


27, 


.693 


0, 


.491 


0, 


.360 


16, 


.831 


14, 


.568 


-4. 


.0 


0, 


.082 


0. 


.071 


0, 


.0 


19, 


.9 


ELAISR142701+331222B 


2.324 


0, 


.018 


2.858 


0, 


.022 


14 


27 


01, 


.5165 


+33 


12 


13, 


.133 


0, 


.502 


0, 


.371 


15, 


.773 


15, 


.202 


86. 


.1 


0, 


.124 


0, 


.119 


0, 


.2 


20, 


.1 


ELAISR142701+331222T 


3.712 


0, 


.025 


7.433 


0, 


.031 


14 


27 


01, 


.8870 


+33 


12 


22, 


.086 


0, 


.486 


0, 


.351 
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Figure 8. SMALL EXAMPLE (This figure is available in 
its entirety in [ittp: // www. ast .cam. ac . uk / ^ciliegi /elais/ papers /[ ) 
: Contour maps for multiple sources in Nl. Contour levels are 
5,6,7,8,10 .... times the local rms noise value (see Table |2| Each 
map is 1' X 1' except for very extended sources for which the 
map size is variable. For these sources the vertical bar in the map 
shows 1 arcmin. 



Fieure 9. As in Figure ^ 


] for N2. This figure 


http: / /www. ast. cam. ac.uk/ 


~cilicgi / elais /papers/ 



and < A<5 0.18". The sources coordinate reported in 
the catalogue (see Table 3) have been corrected for these 
off-sets. 

From the distribution of < Aa > and < A(5 > shown 
in Figure |^ we have estimated the rms calibration errors: 
Ea = 0.486 and es — 0.351". Their small values show that 
our mosaic maps are not affected by relevant geometric dis- 
tortions due to the approximation of a finite portion of the 
spherical sky away from the instrumental zenith with a bi- 
dimensional plane (see Perley 1989, Condon et al. 1998). 
Using the mean off-sets and calibration errors obtained us- 
ing the 37 sources stronger than 10 mjy in common with 
the FIRST survey, we have calculated the positional errors 
of all our sources. To verify that our positional errors are 
realistic also for faint sources, we used the 211 single com- 
pact sources stronger than 1 mJy found to be in common 
with the FIRST survey. The right ascension and declination 



Figure 10. As in Figure [ s] for N3. This figure is available in in 
http: / /www. ast. cam. ac.uk/ ^^ciliegi/clais/papers/ 
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Figure 11. Position errors for strong point sources (37 sources 
in our survey and in the FIRST survey with fiux radio greater 
than 10 mJy) 



differences < Aa > and < AS > were divided by the uncer- 
tainties (Tq and as (from Equations 3 and 4). As expected, 
the normalized error distributions (Figure [l2|) are Gaussian 
with nearly zero mean and unit variance, verifying that our 
positional uncertainties are accurate also for sources down 
to ~1 mJy. 

However, as shown in the flux distribution of Figure ^ 
the majority of our radio sources have a flux lower than 
1 mJy. Since we do not have other radio catalogues below 
1 mJy to repeat out test, we decided to use our observa- 
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Figure 12. Position errors for tiie 211 point sources stronger 
than 1 mjy in common with the FIRST survey, in unit of the 
calculated position uncertainties and cr^, along with the 90% 
confidence error circle [upper panel). The smooth curves in the 
central and lower panel represent the expected Gaussian of zero 
mean and unit variance. 



tions to test the reliability of the position errors down to 
the flux limit of our survey (~ 0.135 mJy). In particular, 
we used the following procedure: (1) Each mosaic map has 
been obtained combining five or more individual pointings, 
whose centers are separated by ~ 22 arcmin. (2) Because in 
each individual map the FWHM of the primary beam is 31 
arcmin, there are many overlapping regions where sources 



are detected in two different pointings. (3) We considered 
the sources detected in each individual pointing as an inde- 
pendent data set and we used the sources in common (i.e 
the sources detected in the overlapping regions) to test the 
reliability of the positional errors. We have a total of 134 
sources in the overlapping regions with a flux lower than 1 
mJy. The right ascension and declination differences of these 
134 sources were divided by the combined errors auA and 
o-DBc, where = cr^i + 0-^2 and ajj^c = ^si + <^^2 with 

CTq2, o"5i and as2 from Equations 3 and 4. As shown 
in Figure |l^, also these error distributions for faint sources 
agree with the expected Gaussians (smoothed curves). 

In conclusion the positional uncertainties obtained us- 
ing Equations 3 and 4 are accurate also for sources down to 
the limit of our survey (~ 0.135 mJy) and they can be used 
to estimate the reliability of the optical and infrared identi- 
fications. The typical rms position uncertainties Oa and as 
are plotted as function of the flux density in Figure |l^. The 
positional errors of the radio sources are ~ 2 arcsec for the 
fainter sources (~0.13 mJy) and ~ 0.6 arcsec for the brighter 
sources (> 10 mJy). 

5.1.2 Bandwidth Smearing 

The principles upon which synthesis imaging are based are 
strictly valid only for monochromatic radiation. When ra- 
diation from a finite bandwidth smearing is accepted, aber- 
rations in the image will result. These take the form of ra- 
dial smearing which worsen with increased distance from the 
center. The peak response to a point sources simultaneously 
declines in a way that keeps the integrated flux constant. 
In other words the bandwidth smearing reduces the peak 
flux density Sp of a source but not its integrated flux den- 
sity S/. However, in our case, since we used the spectral 
line mode where the bandwidth is narrower than in contin- 
uum mode, the bandwidth smearing should be negligible. To 
verify that this is the case, a strong calibrator (the source 
3C286) was observed at position offset of 5,10,15,20,25,30,35 
arcmin in two orthogonal directions (North-South and East- 
West) and we made a Gaussian fit of the source in each of 
the offset position. The ratio between peak and total flux 
density (Sp/S/) as function of the off axis value is shown 
in Figure |l^. As shown in Figure ^ the bandwidth smear- 
ing is negligible up to r~30 arcmin from the center of each 
single pointing. Therefore the resulting mosaic maps, where 
the distance between the single pointings is ~ 22 arcmin, 
are not affected by relevant bandwidth smearing and point 
sources should have Sp/S/ >0.98 everywhere. 



6 THE SOURCE COUNTS 

The complete sample of 867 sources with Speak > 5a was 
used for the construction of the sources counts. Sources with 
multiple components are treated as single radio sources. For 
all the sources we have used the total flux in computing the 
source counts. Every source was weighted for the recipro- 
cal of its visibility area (yi{Sp)~^ , see Figure ^ that is the 
area over which the source could have been seen above the 
adopted limit. The 1.4 GHz source counts of our survey are 
summarized in Table For each flux density bin, the av- 
erage fiux density in each interval, the observed number of 



© 1998 RAS, MNRAS 000, 0-|l| 



12 P. Ciliegi et al. 





6 




4 


o 




LiJ 


2 






b 





'O 


-2 


<1 






-4 




-6 




-4 -2 2 4 6 

A a / CT RA 




■10 -5 5 10 

A a / CT RA 




Figure 13. Position errors for tlie 134 point sources in the over- 
lapping region with flux lower than 1 mjy, in unit of the combined 
position uncertainties cr^^ and Cij^Ci along with the 90% confi- 
dence error circle (upper panel). The smooth curves in the central 
and lower panel represent the expected Gaussian of zero mean 
and unit variance. 



sources in each flux interval, the differential source density 
(in sr~^ Jy~^), the normalized differential counts nS^'^ (in 
sr^^ Jy^'^) with estimated errors (as ri^^'^S'^'^) and the in- 
tegral counts (in sr^^) are given. The normalized 1.4 GHz 
counts are plotted in Figure |l^ where, for comparison, the 
differential source counts obtained with other 1.4 GHz ra- 
dio surveys are also plotted while the integral source counts 
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Figure 14. The typical rms position uncertainties (t„ and cFj for 
point sources of fiux-density Speafe- 
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Figure 15. Bandwidth smearing reduces the peak flux density 
Sp of a point sources but preserves the integrated flux density Sj. 
Abscissa : angular distance from the center (arcmin). Ordinate ; 
The observed ratio Sp/Sj 



(deg~^) are plotted in Figure ^ The solid line in Figure |l^ 
represents the global fit to the counts obtained by Wind- 
horst, Mathis & Neuschaufer (1990) by fitting the counts 
from 24 different 1.4 GHz surveys. The open circles repre- 
sent the counts obtained using the 1.4 GHz survey FIRST 
(White et al 1997). Finally the open stars represent the 1.4 
GHz counts obtained by Gruppioni et al. (1997) in a recent 
radio survey in the Marano Field (RA=03'' 15", DEC=-55 
13). 

As shown in Figure, there is a very good agreement 
between our counts and those obtained with other surveys. 
In particular, the points at the fainter ffux level confirm 
the well-know flattening observed in the normalized differ- 
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ential source counts below few mjy (Windhorst, Mathis 
& Neuschaufer (1990)). The roll-off of the sources in the 
FIRST survey (open circles in Figure ^) at flux densi- 
ties less than ~2 mJy is due to the low peak flux densities 
of many faint, extended sources, which make them unde- 
tectable in the FIRST survey. As noted by White et al. 
(1997) this is a clear indication of the incompleteness of the 
FIRST survey at the faint limit. On the other hand, the 
very good agreement between our counts and those obtained 
with other surveys (also at very faint flux level) and the lack 
of any roll-off at faint flux level in our counts indicate that 
our procedure for the source extraction (see above) yields 
very good results and that our sample is not affected by 
incompleteness at the faint limit. 

A Maximum Likelihood flt to our 1.4 GHz counts with 
two power laws: 



1000 



dN 

Is 



S-^i ffS>S6 
S-°'^ ff S < St 



(5) 



gives the following best fit parameters: qi = 1.71 ± 0.10, 
02 = 2.36 ± 0.20, Sb ~0.5 mJy. These values suggest 
that the re-steepening of the integral counts toward an Eu- 
clidean slope start below 1 mJy, in agreement with Gruppi- 
oni et al. (1997), Condon & Mitchell (1984) and Windhorst, 
van Heerde & Katgert (1984), while Windhorst, Mathis & 
Neuschaufer (1990), by fitting the counts of several 1.4 GHz 
survey, found that the change in slope starts around 5 mJy. 



7 COMPARISON WITH THE NVSS AND 
FIRST RADIO CATALOGUES 

The regions that we are observing with the VLA are covered 
also by the NVSS and FIRST radio survey. The NVSS 
(NRAO VLA Sky Surveys) covers the sky north of J2000 
5=— 40°at 1.4 GHz with 9=45" resolution and nearly uni- 
form sensitivity of 0.45 mJy (1 a), while the FIRST (Faint 
Images of the Radio Sky at Twenty cm) covers over 10,000 
squares degrees at 1.4 GHz with a typical 1 a rms of 0.15 
mJy and a resolution of 5.4". 

It was therefore natural to make a comparison between 
our results and the results obtained by the NVSS and 
FIRST surveys on the same regions of the sky. 



7.1 SAD on the NVSS and FIRST maps 

As first step we used the NVSS and FIRST radio maps 
to test the reliability of the software that we are using to 
extract the radio sources (SAD in the AIPS version October 
96). From the NVSS and FIRST archive we retrieved the 
maps covering the region of the sky in the field N2 covered 
also by the pointing N2 VI, N2 V2 and N2 V3 of our VLA 
observations (see Table ^). 

From the NVSS archive we retrieved the map 
I1640P40 while from the FIRST archive we retrieved 
the maps F16330+40417, F 16330-^4 11 32, F16360+40417, 
F16360+41132, F 16390-1-404 17 and F16390-K41132. 

Using the software SAD on these maps, we obtained two 
lists of sources (one for the NVSS and one for the FIRST) 
in the same region covered by our VLA observations. These 
two lists were compared with the list of sources obtained 
from the NVSS and FIRST catalogs in which are reported 
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Figure 16. The normalized differential source counts. The ab- 
scissa gives the flux density (mJy) and the ordinate gives the 
differential number of sources normalized by S"^'^ (sr~^ Jy+'^'®). 
The solid line represents the global flt to the counts obtained by 
Windhorst, Mathis and Neuschaufer (1990) by fitting the counts 
from 24 different 1.4 GHz surveys. The open circles represent the 
counts obtained using the 1.4 GHz survey FIRST (White et al 
1997), the flUed circles are the counts obtained using our VLA 
observations in the ELAIS regions while the open stars represent 
the 1.4 GHz counts obtained by Gruppioni et al. (1997) in the 
Marano Field. 



the sources above the 5cr limits of each survey. The results 
are reported in Table ^ and in Figure 

All the three NVSS sources missing in the list obtained 
by us with SAD have a peak flux density (calculated by us 
with MAXFIT on the NVSS map) lower than ~ 2.0 mJy (i.e 
< the 5 a limit of the NVSS maps). Therefore they are 
probably spurious sources. 

Three (163637.9+410511, 163707.54-405125 

and 163815.4-1-405840) of the eight sources present in the 
FIRST catalogue but not detected by us using SAD on the 
FIRST map are not detected also in our deeper VLA maps 
to a flux limit of 0.1 mJy, so they are spurious sources. The 
other five sources missing in the FIRST list (as well as the 
four sources detected by SAD but not in the FIRST cata- 
log) have a peak fiux density < 1.2 mJy, i.e. very close to 
the FIRST fiux limit. These differences are probably due 
to the unreliable Gaussian fit of faint sources as discussed 
above. Moreover we have to consider that the sources lists 
(our SAD list and the catalogue lists) are obtained using dif- 
ferent source extraction software. In fact, both NVSS and 
FIRST sources are extracted with special AlPS-based rou- 
tines wrote exclusively for each survey (see Condon et al. 
1993 and White et al. 1997 for more details). 

Therefore, in conclusion, we are confident in using SAD 
as source extraction software, keeping in mind that for faint 
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Table 5. The 1.4 GHz Radio Source Counts 
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53 


1.166x10'^ 
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59.07± 14.33 
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13 
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162. 10± 57.31 


1.166x10-' 


83.55 


- 150.39 
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Table 6. Comparison between SAD and NVSS - FIRST catalogs 



SURVEY 


Catalogue 


SAD 


Common 


Catalogue 


SAD Sources 




Sources 


Sources 


Sources 


Sources 


not in Catalogue 










not in SAD 




NVSS 


46 


43 


43 


3 





FIRST 


74 


70 


66 


8 


4 



sources the flux obtained with a Gaussian fit may be unre- 
liable. 



Table 7. Sensitivity Comparison NVSS, FIRST and our VLA 

survey 



7.2 Comparison between our VLA observations 
and the NVSS and FIRST source catalogues 

The primary difficult in comparing the three surveys is the 
difference in angular resolution. A summary of the main 
properties of the NVSS, FIRST and our VLA survey is re- 
ported in Table ^. 

The high resolution surveys will partially resolve ex- 
tended sources that appear point-like in the low resolution 
surveys, so that the flux density of the sources will be smaller 
in FIRST and our survey than in the NVSS. To mini- 
mize this problem, we restricted the comparison to compact 
sources {i.e sources with S//Sp <2) with an off-axis value 
in our VLA map lower than 34 arcmin. The restriction on 
the off-axis value allow us to use only the sources detected 
in regions where the 5a limit is well below the limit of the 
FIRST survey (see Table | and Table 0). For the NVSS 
and FIRST surveys we used the public source catalogues 
as available on February 1, 1998. 

7.2.1 Our VLA survey versus the NVSS survey 

Using the restriction of an off-axis value lower than 34 ar- 
cmin and a peak flux density Sp >2.3 mJy (~ the 5<j rms 





NVSS 


FIRST 


Our VLA 

data 


VLA Configuration 


D 


B 


C 


Beam size 


45" 


5" 


15" 


confusion noise 


0.08 mJy 


0.002 mJy 


0.011 mJy 


Exposure time 


30-60 s 


3 min 


~1 hour 


la rms noise 


0.45 mJy 


0.15 mJy 


0.05 mJy 


5cr rms limit 


2.25 mJy 


0.75 mJy 


0.25 mJy 



Note: The rms confusion from faint extragalactic source in 
C-configuration at 20cm is O.OlmJy (Mitchell & Condon, 1985). 
The value for D-configuration comes from Condon et al. 1998. 

limit of the NVSS survey) we have 109 compact sources in 
common with the NVSS survey. In Figure ^ we show a com- 
parison between our VLA and NVSS total flux densities. A 
comparison between the flux scales is discussed below. Be- 
sides these 109 common sources, there is one NVSS source 
(NVSS 142753+330837) that we did not detected in our sur- 
vey plus 5 sources (with Sp >2.3 mJy) detected in our maps 
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Figure 17. The integral source counts. The abscissa gives the 
flux density (mJy) and the ordinate gives the total number of 
sources per square degree 

but not in the NVSS survey. In Figure ^ and ^ we report a 
contour plot of these sources. For comparison a contour plot 
of the FIRST image in the same region of the sky is also 
reported. As shown in Figure ^the NVSS sources missing 
in our survey is not detected also in the FIRST survey. Its 
NVSS radio flux is 2.3 mJy, very close to the 5 a detection 
limit of the survey. It is probably a spurious source but it 
could be also a very low surface brighteness source detected 
only in the NVSS, the radio survey with the lowest reso- 
lution. However a strange phenomenon like a variable radio 
source can not be completely excluded. For the 5 sources 
missing in the NVSS survey, Figure ^ clearly shows that 
one (ELAISR142917-I-332626) is missing because our higher 
resolution survey has resolved an extended NVSS source, 
one (ELAISR142940+330552) is completely undetected on 
the NVSS map while the other 3 sources are detected also 
in the NVSS map but below the 5 a limit (~ 2.25 mJy). 




10 100 

NVSS Catalog Flux DcnsiLy (rnjy) 



FIRST sources in N2 



FIRST Catalog Mux Density (iTiJy) 

Figure 18. Ratio of flux reported in the public catalogue to our 
flux obtained using SAD on the radio map as function of the flux 
reported in the catalogue for the NVSS (top panel) and for the 
FIRST (lower panel) surveys 



All these five sources are present in the FIRST cat- 
alogue (and in the FIRST map as clearly showed in Fig- 
ure ^l|). Their absence in the NVSS catalogue is an indi- 
cation of the incompleteness of the NVSS survey near the 
flux limit. 



7.2.2 Our VLA survey versus the FIRST survey 

Using the restriction of an off-axis value lower than 34 ar- 
cmin and a peak fiux density Sp >1.0 (the limit of the 
FIRST catalogue), we have 215 compact sources in common 
with the FIRST survey. In Figure we show a comparison 
between our VLA and FIRST total flux densities. 

Besides these 215 common sources, there are 14 FIRST 
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Figure 21. Contour plot of the 5 sources detected in our VLA and FIRST survey but not in tlie NVSS survey. All the maps have a 
size of 3' X 3' and a the contour levels of 1.5,1.8,2.0,2.2,2.5,3,5,8 mjy. 
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Figure 19. The flux obtained with our VLA survey vs. the 
NVSS flux density as reported in the pubUc catalogue (upper 
panel) and the ratio of our VLA to NVSS flux density as func- 
tion of the NVSS flux density (lower panel). 



sources that we did not detected in our survey and 29 
sources present in our survey but not in the FIRST cat- 
alogue. A flux distribution and a contour maps of these 
last 29 sources are shown in Figure ^ and |24| while 
a contour maps of the 14 FIRST sources that we did 
not detect are shown in Figure ^ As shown in Fig- 
ure ^ all the radio sources missing in the FIRST cat- 
alogue have a peak flux density lower than 2 mJy. How- 
ever, as shown in Figure many of them appear on the 
FIRST maps. This result confirms the incompleteness of 
the FIRST survey below ~2 mJy, as already indicated by 
the differential source counts reported in Figure ^. On the 
other hand, Figure ^ shows that many of the 14 FIRST 
sources missing in our survey probably are not real. Some of 
them may be uncleaned residual around strong sources (see, 
for example, FIRST161318-(- 541607, FIRST163646+405442, 
FIRST163815+405840) or simply spurious sources in a 
not well cleaned map (see, for example, the FIRST 
maps of FIRST161351-H543258, FIRST163634-h413013 and 
FIRST163706-H405125). 



7.2.3 Flux comparison between the three surveys 

As shown in Figure |l^ and ^ the flux densities of our 
survey are in good agreement with the FIRST and NVSS 
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Figure 20. Radio contour plot of the NVSS source NVSS 
142753 +330837 present in the NVSS catalogue but not in our 
VLA and FIRST catalogue overimposed to the optical Digi- 
tized Sky Survey (DSS) image. All the maps have a size of 
3' X 3'. Radio contour level are 0.25,0.35,0.5,0.75,1.0,1.5,1.8,2.5 
mJy for our map, 0.75,1.0,1.5,1.8,2.5 mJy for the FIRST maps 
and 1.0,1.3,1.5,1.7,1.8,2.0,2.2,2.5 mJy for the NVSS map 



flux densities over two order of magnitude. However, as ex- 
pected, the high resolution surveys tend to estimate lower 
flux than the lower resolution survey. This effect is evident 
in the lower panels of Figure ^ and ^ : our VLA flux 
densities are, in mean, lower than the NVSS flux densities 
but higher than the FIRST flux densities. High resolution 
surveys, with their smaller synthesized beam size, lose flux 
due to the resolution surface brightness effect. However some 
consideration cam be made from Figure |l| and m In spite 
a factor of 3 in angular resolution between our VLA and 
NVSS survey (15" vs. 45" FWHM) the flux ratio of the 
two surveys is always lower than 0.5 while the flux ratio be- 
tween our VLA and FIRST surveys reaches values of ~ 4, 
although the two surveys have still a factor of 3 in angu- 
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Figure 22. The flux obtained with our VLA survey vs. the 
FIRST flux density as reported in the pubhc catalogue (upper 
panel) and the ratio of our VLA to FIRST flux density as func- 
tion of the FIRST flux density (lower panel). 
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Figure 23. Flux distribution of the 29 sources detected in our 
survey with a peak flux density Sp > 1 mJy but not present in 
the FIRST catalogue 



lar resolution differences (5" vs. 15" FWHM). The missing 
flux between the FIRST (VLA - B conflguration) and our 
survey (VLA - C configuration) is greater than the missing 
flux between our survey and the NVSS survey (VLA - D 
configuration). Therefore the C conflguration of our survey 
with a synthesized beam size of 15" seems to be the bet- 
ter compromise between high (B configuration) and low (D 
configuration) resolution radio surveys. It is less prone to 
surface brightness efi^ects than the B configuration without 
an excessive loss of flux in comparison with the D configu- 
ration. 



8 SUMMARY 

Using the Very Large Array (VLA) radio telescope, we 
observed at 1.4 GHz a total area of 4.222 deg^ in the 
ISO/ELAIS regions Nl N2 and N3. The lower flux density 
limit reached by our observation is 0.135 mJy (at 5 a level) 
on an area of 0.118 deg^, while the bulk of the observed re- 
gions are mapped with a flux density limit of 0.250 mJy (5 
a). The data were analyzed using the NRAO AIPS reduction 
package. The source extraction has been carried out with the 
AIPS task SAD. The reliability of SAD has been tested using 
the maps of the radio surveys FIRST and NVSS. 

Considering all the available observations, we detected 
a total of 867 sources at 5 cr level, 44 of which have mul- 
tiple components. These sources were used to calculate the 
normalized differential source counts. They provide a check 
on catalogue completeness and reliability plus information 
about source evolution. A comparison with other surveys 
shows a very good agreement, conflrming the presence of 
the well-know flattening of the counts below 1 mJy, the com- 
pleteness of our catalogue and the reliability of our proce- 
dure for the source extraction. 

A comparison with the FIRST and NVSS radio sur- 
veys has conflrmed the incompleteness of these two surveys 
near their flux limits, while a flux comparison between the 
three surveys has shown that our survey with the VLA ar- 
ray in C conflguration is the best compromise between high 
and low resolution radio surveys. The positional errors of the 
radio sources are ~ 2 arcsec for the fainter sources (~0.13 
mJy) and ~ 0.6 arcsec for the brighter sources (> 10 mJy). 
This small value will enable us to obtain an accurate and 
fast optical/infrared identiflcation of the radio sources. 
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Figure 24. SMALL EXAMPLE (This figure is available in 
its entirety in [ittp: / / www. ast .cam. ac . uk / ^ciliegi /elais/ papers /| ) 
: Contour plot of the sources detected in our VLA survey with 
a peak flux density Sp > 1 mjy but not present in the FIRST 
catalogue. The dashed circle shows the position of the source de- 
tected in our VLA survey. All the maps have a size of 2' X 2' and 
a contour levels of 0.50,0.8,1,1.5,2,5,8,10 mJy 



Figure 24. SMALL EXAMPLE (This figure is available n 
its entirety in http:/ /www. ast. cam. ac.uk/~ciliegi/elais/papers/): 
Contour plot of the FIRST source not in our VLA catalogue. 
The dashed circle shows the position of the FIRST source. 
All the maps have a size of 2' X 2' and a contour levels of 
0.50,0.8,1,1.5,2,5,8,10 mJy. 
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